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High-resolution spectroscopy for Cepheids distance determination 

II. A period-projection factor relation 
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ABSTRACT 

Context. The projection factor is a key quantity for the interferometric Baade-Wesselink (hereafter IBW) and surface- 
brightness (hereafter SB) methods of determining the distance of Cepheids. Indeed, it allows a consistent combination 
of angular and linear diameters of the star. 

Aims. We aim to determine consistent projection factors that include the dynamical structure of the Cepheids' atmo- 
sphere. 

Methods. Hydrodynamical models of 5 Cep and I Car have been used to validate a spectroscopic method of determining 
the projection factor. This method, based on the amplitude of the radial velocity curve, is applied to eight stars ob- 
served with the HARPS spectrometer. The projection factor is divided into three sub-concepts : (1) a geometrical effect, 
(2) the velocity gradient within the atmosphere, and (3) the relative motion of the "optical" pulsating photosphere 
compared to the corresponding mass elements (hereafter /o-g). Both, (1) and (3) are deduced from geometrical and 
hydrodynamical models, respectively, while (2) is derived directly from observations. 

Results. The Fe I 4896.439 A line is found to be the best one to use in the context of IBW and SB methods. A co- 
herent and consistent period-projection factor relation (hereafter Pp relation) is derived for this specific spectral line: 
p = [—0.064 ± 0.020] logP + [1.376 ± 0.023]. This procedure is then extended to derive dynamic projection factors for 
any spectral line of any Cepheid. 

Conclusions. This Pp relation is an important tool for removing bias in the calibration of the period-luminosity relation 
of Cepheids. Moreover, it reveals a new physical quantity /o-g to investigate in the near future. 

Key words. Techniques: spectroscopic - Stars: atmospheres - Stars: oscillations (including pulsations) - (Stars: vari- 
ables): Cepheids - Stars: distances 



K^ 1. Introduction In 1993, Butler studied the velocity gradient in the at- 

k^ mosphere of four Cepheids using the excitation potential of 

'V^ ' spectral lines, together with their asymmetry. Then, Butler 

^ , The period-luminosity relation (hereafter PL relation) of g^ ^1. (1996) introduced the velocity gradient in hydrostatic 

■ - - ■ the Cepheids is the basis of the extragalactic distance scale, g^gHar atmosphere models and found a 20% reduction on 

but its calibration is stiU uncertain at a AM = ±0.10 mag ^^^ amplitude of the pulsation velocity curve in the case of 

level. Long-baseline interferometers currently provide a ^ ^ql. In addition, the 7-velocity was reduced by 2 km s-\ 

new, quasi-geometric way to calibrate the Cepheids PL re- „ , . ^ , /,„„„n _l i- 1 _li 1 -^ i- ^ • ^i 

1 , . T J ° , ., . -ui 4- J 4. • 4-1, J- 4- l-okm et al. (1996) studied the velocity gradient m the 

lation. Indeed, it is now possible to determine the distance , , ,. „^„ ' , , , , ■ , , ,t 

r 1 ,• r~i 1 -J i 11 -4.1 4.1 • J. r 4. • atmosphere ot Cep, based on nydrodynamical modelling. 

01 galactic Cepheids up to Ikpc with the mterterometric ^^ ^^ ^ ^ ^ ■ ^ f r 1 1 • r 

r, 1 ^TiT 1-1 ii 1 1 Ti TT1TT7 4_i 1 n It was touud to DC au important source oi broadening tor 

Baade-Wessehnk method, hereafter IBW method; see tor ^ ... ^ , ,. , .^ . ., ^ ^ ^. *= ^ 

e.g. Sasselov & Karovska (1994) and Kervella et al. (2003). '^''^f''' ^"^l^^'^^ J^^f ' ^^^^^ ^^^l^"" *° "^"^^^^^^ ">' g^°™^*- 

T° r , . ^ 1 J 4- 1 J- 4. rical proiection enects. 

Intcrteromctric measurements lead to angular diameter es- 

timations over the whole pulsation period, while the stellar Using an improved version of the hydrodynamical model 
radius variations can be deduced from the integration of the of S Cep, Nardetto at al. (2004, hereafter Paper I) pro- 
pulsation velocity. The latter is linked to the observational POsed an interferometric definition of the projection factor, 
velocity deduced from spectral line profiles by the projec- ^ difference of a few km s ^was found between the ampli- 
tion factor p. In this method, angular and linear diameters tude of the photospheric and line-forming region's velocity 
have to correspond to the same physical layer in the star curves, leading to a bias of 5% on the derived distance. This 
to correctly estimate the distance. The projection factor is theoretical result has been observationally confirmed using 
currently the most important limiting quantity of the IBW trigonometric parallax measurements of the HST and opti- 
method. Indeed, in addition to limb-darkening effects, it is cal long basehne interferometry by Merand et al. (2005). 
related to the velocity gradient and, more generally, to the Nardetto et al. ([20063) show that spectro- 
dynamical structure of the Cepheid atmosphere. interferometry provides a new geometric view of the 
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Cepheids' atmosphere. However, the combination of differ- 
ent techniques (high-resolution spectroscopy, spectro- and 
differential- interferometry) is needed to efhciently con- 
strain the physical parameters of the Cepheid atmosphere 
and, in particular, the projection factor. 

Recently, while comparing radial velocity curves of dif- 
ferent species, Petterson et al. (2005) found some evidence 
of a relation between the velocity gradient in the Cepheids' 
atmosphere and their period. Using a selected sample of 
absorption metallic lines, we propose to probe the veloc- 
ity gradient in the Cepheids' atmosphere in order to de- 
termine dynamic projection factors. First, by using 5 Cep 
and I Car hydrodynamic models, we present a new spectro- 
scopic method for determining the velocity gradient. Then, 
the method is applied to the eight Cepheids observed with 
the HARPS instrument. We discuss the choice of the spec- 
tral line and then derive a specific period-projection fac- 
tor relation (hereafter Pp relation). Results are discussed 
within the framework of the PL relation. Finally, we pro- 
pose a general method (for all lines and all stars) of deter- 
mining the projection factor. 



Table 1. Spectral lines used in this study. 



Name 


Wavelength (A) 


Fel 


4683.560 


Fel 


4896.439 


Fel 


5054.643 


Nil 


5082.339 


Fel 


5367.467 


Fel 


5373.709 


Fel 


5383.369 


Till 


5418.751 


Fel 


5576.089 


Fel 


5862.353 


Fel 


6024.058 


Fel 


6027.051 


Fel 


6056.005 


Si I 


6155.134 


Fel 


6252.555 


Fel 


6265.134 


Fel 


6336.824 



2. Cepheids observed and selected spectral lines. 

Ten stars have been observed with the HARPS spectrom- 
eter {R = 120000) : R Tra, S Cru, Y Sgr, (3 Dor, C Gem, 

Y Oph, RZ Vel, I Car, RS Pup, and X Sgr. In the first 
paper of this series, Nardetto et al. (2006b, hereafter paper 
H) showed that the radial velocity associated with the cen- 
troid of the spectral line, together with the line asymmetry, 
are very important tracers of the dynamical structure of 
the Cepheids' atmosphere. X Sgr was studied separately by 
Mathias et al. (2006) because of its very atypical behavior 
showing several components in the spectral lines profiles. 

Y Oph is not considered here due to its insufficient phase 
coverage (see Paper H, Fig. 3). 

Using Kurucz's models ()1992[ ). we identified about 150 
unblended spectral lines. In Paper II, we considered only 
the unblended metallic line Fel 6056. 005A. In this second 
paper, we have carefully selected 17 spectral lines following 
two criteria: (1) in order to avoid bias in the determination 
of the line depth, the continuum must be perfectly defined 
for all pulsation phases and for all stars. An example of the 
quality required is given in paper II for the Fel 6056.005 
A spectral line (see Fig. 1); (2) the selected sample of lines 
has to cover a wide range of depth. The selected spectral 
lines are presented in Table[TJ Depending on the star and 
the spectral line considered, the line depth can range from 
2% to 55%. 

3. Hydrodynamical models. 

The hydrodynamical model of 5 Cep is presented in Paper I. 
In addition, we have derived a new, consistent model of 
£ Car. Since the main stellar quantities of I Car (HD 84810) 
are still uncertain, we tried several sets of luminosity L, 
effective temperature Tcft, and mass M in order to get suit- 
able observational quantities, such as the pulsation period 
P, the average radius of the star i?, bolometric and radial 
velocity curves, and the line profiles. The ML relation was 
taken from Chiosi et al. (1993), and the OPAL opacity ta- 
bles (Rogers & Iglesias (1992)) were used. 

This leads to the following set of parameters for a 154- 
zone model: M = 11.5 M©, L = 21000 L©, T^e = 5225i^, 



Table 2. Hydrodynamical models of Cepheids. 



Name 


P 

[days] 


Teff 

[K] 


L 

^0 


M 

Mq 


R 

Rq 


S Cru 


4.7 


5900 


1900 


5.6 


42 


<5Cep 


5.4 


5877 


1995 


4.8 


43 


YSgr 


5.7 


5850 


2200 


5.0 


45 


/3Dor 


9.9 


5500 


3500 


5.5 


65 


C Gem 


10.4 


5500 


3600 


5.0 


64 


RZ Vel 


21.6 


5400 


7450 


7.0 


109 


f Car 


34.4 


5225 


21000 


11.5 


180 


RSPup 


42.9 


5100 


22700 


9.7 


186 



Y = 0.28, and Z = 0.02, which corresponds to a typical 
Pop. I chemical composition. The inner boundary has been 
fixed at about T — 3.5 • 10^ K, corresponding to 4% of 
the photospheric radius, so the model envelope with the 
atmosphere contains about 53% of the stellar mass. The 
atmosphere itself contains about 10^'' of the total stellar 
mass. 

We started the hydrodynamical calculations with a lin- 
ear, non-adiabatic fundamental-mode velocity profile hav- 
ing a value of 10 km s~^ at the surface. At the limit cycle, 
the pulsation period is 34.4 days, very close (3%) to the 
observational (P = 35.6 days) value deduced by Szabados 
et al. (|l989f ). Bolometric and radial velocity amplitudes are 
1.3 mag and 50 km s^^, respectively. The relative radius 
amplitude at the surface is AR/R — 17%. The mean pho- 
tospheric radius is about R — 180 Pq. 

For the hydrodynamical models of S Cep and £ Car used 
in this paper (see Sect. |4]), a careful analysis of the dynam- 
ical structure of their atmospheres was performed based on 
radiative transfer computations (under local thermal equi- 
librium) of all spectral lines in Table[T] 

We only derived the photospheric pulsation velocity for 
other Cepheids (see Sect. 14. 4| ). Thus, the physical parame- 
ters are only roughly estimated since the dynamical struc- 
ture of the Cepheids' atmosphere is not considered. The 
period and the radius of the Cepheids are consistent with 
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observations at the 2% and 3% levels respectively, while 
the mass follows the period-mass relation of Choisi et al. 
(1993) within 13%. The effective temperature and luminos- 
ity were mostly constrained by the period and the radius. 
R TrA was not modeled because of its extremely short pe- 
riod. A specific and in-depth study would be necessary to 
model this star. 

The hydrodynamical models are presented in Table [2l 



4. New insights into understanding the projection 
factor 

In this section, we propose a division of the projection factor 
into sub-concepts in order to allow a direct constraint from 
HARPS spectroscopic observations. To test this method, 
we consider the hydrodynamic models of S Cep and £ Car. 



4.1. The projection factor definition 

First of all, we have to provide a definition for the projection 
factor that should be applied in the IBW and SB methods. 
This has already be done in Paper I, but it must now be 
refined and adapted to the method proposed here. 
We define the interferometric projection factor as 



AVl_ 



(1) 



where AVp° is the amplitude of the pulsation velocity curve 
associated to the photosphere (subscript p) of the star. It 
corresponds exactly to the optical (subscript o) barycenter 
of the photosphere defined by Tc — 2/3, where t^ is the op- 
tical depth in the continuum. ARVc is the amplitude of the 
radial velocity curve obtained with the centroid method, 
i.e. the first moment of the spectral line. This definition is 
justified for the following reasons: 

First, we consider the pulsation velocity instead of the 
radius, as already proposed in paper I, in order to allow a 
direct application to spectroscopic observations. 

Second, we consider velocity amplitudes to avoid dif- 
ficulties related to the 7-velocity. The 7-velocity is the 
averaged value of the radial velocity curve over one pul- 
sation period. This quantity depends on the line consid- 
ered. Moreover, as shown in Paper I, the projection factor 
is mainly constrained by velocity amplitudes, which also 
justifies this choice. 

Third, we consider the RVc velocity instead of the radial 
velocity associated to the Gaussian fit method, as in Paper 
I. This is required for obtaining a projection factor inde- 
pendent of the rotation of the star and the natural width 
of the spectral lines (see Fig. 8, Paper II). We insist on 
this definition of the radial velocity since it is absolutely 
required to allow important comparisons between the pro- 
jection factors of Cepheids. 

In the next section, we consider the Fel 4896.439 A 
spectral line as a reference. But the results, in terms of 
consistency, can be generalized to any other spectral line. 
We find p = 1.33 for S Cep and p = 1.27 for £ Car. These 
values are our references in the following. If we apply a min- 
imization process between V°{4>) and pRVc{(f>), with p the 
only free parameter, we find differences of 0.01, which pro- 
vides a good estimate of the uncertainty on the projection 
factor. 



4.2. Decomposition of tlie projection factor 
We now divide the projection factor: 

P^Po /grad /o-g (2) 

into different quantities where, 

— /o-g = x^, where AV^ is the gas (subscript g) ve- 
locity associated to the optical barycenter (tc = 2/3) of 
the photosphere. Thus, /o-g is linked to the distinction 
between the optical and ga.s photosphcric layers. The 
optical layer is the location where the continuum and 
line photons are generated (e.g. the location of the pho- 
tosphere) . The gas layer is the location of some mass el- 
ement in the hydrodynamic model mesh where, at some 
moment in time, the photosphere is located. Given that 
the location of the photosphere moves through differ- 
ent mass elements as the star pulsates, the two "layers" 
have different velocities, hence the necessity of the /o-g 
definition. Indeed, the interferometer in the continuum 
is only sensitive to the optical layer. 

^ /grad — Av^' whcrc AT^j^ is the gas velocity associated 

to the optical barycenter {t\ = 2/3) of the line-forming 
(subscript I) region. Thus, /grad is linked to the velocity 
gradient in the atmosphere of the star. This quantity 
depends on the line considered. 

— Po = A ny is the geometrical projection factor. It corre- 
sponds to an integration of the pulsation velocity field 
(associated with the line-forming region) projected on 
the line of sight and weighted by the surface bright- 
ness of the star (including limb-darkening in the spec- 
tral line). To derive Po, we use intensity distributions 
in the continuum p rovided by the model. As described 
in Nardetto et al. (|2006^, there is a relation between 
the pulsation velocity and the limb-darkening (and thus 
Po): the intensity profile corresponding to the highest 
velocity at contraction is the most limb-darkened, while 
the profile corresponding to the highest velocity at ex- 
pansion is the least limb-darkened. Given that our def- 
inition of Po is related to the amplitude of radial and 
pulsation velocity curves, we consider the median value 
(peak-to-peak) of thepo(0) curve to derive thepo-factor. 
However, considering intensity distribution in the con- 
tinuum is an approximation. Indeed, hydrodynamic ef- 
fects can result in much limb-darkening, especially at 
the wavelengths corresponding to spectral lines (see e.g. 
Marengo et al. 2003). Nevertheless, this seems to be 
negligible since we obtain a good decomposition of the 
projection factor. 

To test this decomposition, we deduced these three 
quantities directly from the hydrodynamical model, consid- 
ering the Fel 4896.439 A metahic hne. For S Cep, we find 
/o-g = 0.963, /grad = 0.993, and po = 1.390; and for £ Car 
/o-g = 0.944, /grad = 0.982, and Po = 1.366. The corre- 
sponding projection factors are p[S Cep]= 0.963 * 0.993 * 
1.390 = 1.33 and p[£ Car]= 0.944 * 0.982 * 1.366 = 1.27, 
which correspond to our reference values. The physical de- 
composition of the projection factor is thus consistent. 

One can notice that the p-factors presented in this paper 
are independent of the pulsation phase. Indeed, they corre- 
spond to a specific definition based on velocity amplitudes 
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Fig. 1. The projection factor decomposition (p — Po/grad/o-g) in the case of the Fel 4896.439 A spectral line: ARVc (a), 
AFj^ (<), AV^ (v) and AV° (n). Reference values p — 1.33 and p — 1.27 are indicated for 6 Cep and £ Car, respectively. 
The proposed method of the /grad determination is to first derive ARVc as a function of the line depth for all spectral 
lines (lower part, black points). Then, when translating (£), ARVc) points into (I?, PoARVc), a new linear relation is 
found (upper part of the figure, black points). This new relation is (1) coherent with the pulsation velocity gradient in 
the atmosphere (open circles) and (2) its zero-point is consistent with the pulsation velocity corresponding to the gaseous 
layer of the photosphere: AV^ (v). 
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(Eq. [T]). Figure[T] represents the projection factor decompo- 
sition in the case of the Fe I 4896.439 A line for S Cep and 
i Car: ARV^ (a), AV^^ (<), AV^ (v), and AV° (n). 

4.3. Is there a consistent way to derive /giad directly from 
observations? 

Let assume a linear relation between the line depth and 
the position of the line-forming region in the atmosphere of 
the star. Then, considering different spectral lines spread 
all over the atmosphere, it should be possible to determine 
the velocity gradient within the atmosphere. Moreover, by 
extrapolation to the zero line depth, it should be possible 
to reach the photospheric pulsation velocity. 

To test these ideas, we consider the spectral line depth 
corresponding to the minimum extension of the star {D 
in the following). We discuss this choice below. Then, we 
derive ARVc for all spectral lines in Table [1] From this 
quantity, one should be able to probe the velocity gradient 
in the atmosphere and directly derive /grad- For that, we 
plot ARVc as a function of the line depth for all spectral 
lines (Fig.[T|). We find a linear correlation {ARVc — ao-D + 
bo) given by the following relations : ARVc[S Cep]= 2.90I?-|- 
32.84 and ARVc[e Car]= 5.201? -I- 37.02. Interestingly, these 
relations are spectroscopic observables. 

Then, following our decomposition, we translate {D, 
ARVc) points into {D, PoARVc) with po[S Cep]= 1.390 
and Po[i Car]= 1.366 (Fig. [T|). The new relations are then 
: PoARVc[S Cep] = 4.031) + 45.64 and PoARVdl Car] = 
7.10D-I-50.53. Another possible way to derive Po, instead of 
using hydrodynamical modeling, is to consider a linear law 
for the continuum-intensity distribution of the star defined 
by J(cos(0)) = 1 — uv + uvcos(6'), where uy is the limb- 
darkening of the star in the V band (Claret et al. 2000), and 
uy is related to the effective temperature Tgff, the surface 
gravity log (7, the star metallicity, and the turbulent veloc- 
ity. Using this method to derive Po, we find that a slight 
correction must be applied to allow a comparison between 
geometric and hydrodynamic modeling. We find po [hydro] = 
Po[geo]— 0.010 for S Cep and Po[hydro]= Po[geo]— 0.025 for 
i Car. As already mentioned, the limb-darkening in the con- 
tinuum used to derive Po [hydro] is linked to the pulsation 
velocity and, more generally, to the whole dynamical struc- 
ture of the Cepheid atmosphere. Such a dynamical effect 
could explain the difference that we find between po [hydro] 
and Po[geo]. From this, we estimate a 0.01 uncertainty on 
Po- The limb-darkening (and thus Po) will be studied in 
detail by interferometry in the near future. 

We then overplot the amplitude of the modeled pul- 
sation velocity curves corresponding to the different line- 
forming regions AV^^^. Two important points have to be 
mentioned here. 

1. The superimposition oipoARVc and AV^^ values is very 
satisfactory for both stars. Thus, even if the ARVc quan- 
tity includes all the dynamical structure of the line- 
forming region (including the limb darkening), it seems, 
on average, to give direct access to the gas velocity cor- 
responding to Ti = 2/3, i.e. AV^^. In particular, the limb 
darkening (and thus the Po-i&ctor) seems to be indepen- 
dent of the spectral line considered. For any line depth 
D, we stiff have AV,^ = p^ARVc. 

2. We find that the zero-points of the {D, poARVc) rela- 
tions, 45.64 [6 Cep] and 50.53 [£ Car], are close to the ref- 



erence values derived directly from the pulsation veloci- 
ties of the model, AV^ = 45.60[5 Cep] and 50.40[^ Car]. 
This means that, in both cases, the extrapolation to the 
photosphere is verified. We thus have AV^ = Pobo- This 
is the key point that shows that the method proposed 
here is consistent and can be applied to any Cepheid. 
To allow this very important condition, one has to use 
the spectral line depth corresponding to the minimum 
extension of the star {D = D{(j) = i?min))- If another 
estimator is used, for example D = D{4> = Rmax), or 
the line depth averaged over the entire pulsation cy- 
cle D =< D >, then the condition AV^ — Pobo dis- 
appears. Indeed, by extrapolating the three PoARVc 
linear relations corresponding to the three estimators 
{D{(j) = i?inax), <D>, and D{<j) = Emin), respectively) 
towards the zero line depth, we find different agreements 
compared to the reference values derived directly from 
the variation in the photospheric pulsation velocity of 
the model: -1.13%, -0.70%, and 0.09% for S Cep and 
-4.53%, -3.29%, and 0.25% for i Car (see Fig.©. The 
best extrapolation is thus obtained for the 'minimum' 
estimator. We emphasize that this result is verified for 
both models {S Cep and i Car) that have very differ- 
ent physical properties : in period (5.4d, 35. 5d), radius 
(43.4 i?o, 180 i?o), mass (4.8 M©, 11.5 M©), and effec- 
tive temperature (5877-?^, 5225K). The 'minimum' es- 
timator of the line depth (D) is used for all Cepheids in 
the whole paper. Considering this estimator (link to the 
minimum extension of the star) does not mean that our 
projection factors are related to any specific pulsation 
phase. It only means that we have to use this estimator 
to get the right extrapolation of the amplitude (peak-to- 
peak) of the pulsation velocity toward the photosphere. 

Thus, taking the linear correlation ARVc = clqD + 60 
into account, we propose the following relation to observa- 
tionally derive /grad : 



Jgrad 



bo 



(3) 



aoD + bo 

This quantity can be derived directly from observations. 
In the case of our models and using the line depth corre- 
sponding to Fe I 4896.439 A , we find /grad = 0.990 [6 Cep] 
and /grad = 0.988 [£ Car]. These results are very consistent 
with the same quantities derived directly from our projec- 
tion factor decomposition. Consequently, translating these 
quantities into the projection factor p = Po/grad/o-g, we 
find the reference values. We have thus found a new way to 
determine /grad directly from observations. 

4.4. Hydrodynamic modeling of fo-g 

The third quantity /o-g is difficult to determine directly 
from observations. It requires a precise knowledge of the 
dynamical structure of the Cepheid atmosphere and line- 
forming regions. In the cases of 6 Cep and £ Car, the model 
gives /o~g = 0.963 and 0.944, respectively. 

In order to test the Pfo-g relation, we modeled the 
other Cepheids of our sample (see Table [2]) . We only con- 
sider the velocity curves corresponding to the optical and 
gas photospheric layers. Consequently, no radiative transfer 
is calculated, and the resulting physical parameters of these 
models should be considered with caution. Nevertheless, 
these uncertainties are not critical for deriving /o~g. 
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Fig. 2. AVj^ (open circles) as a function of the line depth estimators D(0 = -Rmax), < D >, and D(0 = i?min)- Black 
points correspond to the Po^RVc quantity. The symbol V indicates AVC^; i.e., the amplitude of the photospheric pulsation 
velocity associated to the gas. We find good agreement between AV[^ and PoARVc, which both present an interesting 
linear relation. Through extrapolation (dashed line), we find that the best estimator of the line- forming region is the line 
depth associated with the minimum radius, a result confirmed for 5 Cep and £ Car hydrodynamic models. 



The results are 0.966[S Cru], 0.962[Y Sgr], 0.955[/3 Dor], 
0.953[C Gem], 0.951[RZ Vel], and 0.943[RS Pup]. We con- 
sider a 0.05 uncertainty on fo-g- 

Using these results (based on eight models), we obtain 
the following linear relation: 



/„_g = [-0.023 ± 0.005] log P + [0.979 ± 0.005] (4) 

Using this relation we find fo-g = 0.967 for R TrA. /o-g is 
plotted as a function of the period in Fig. [31 

This linearity can be understood with the following pic- 
ture. Let us assume two atmospheres of short- and long- 
period Cepheids (noted S and L) at expansion. The gas 
velocity is assumed to be the same for both atmospheres. 
Then, due to a geometrical effect (the radius of L is larger 
than the radius of S), the volume of L increases faster than 
S. Consequently, the density and temperature (in the adia- 
batic limit) also decrease faster in L. The lower atomic level 
of this line depopulates faster for L than for S, the opacity in 
the line also falls faster, the spectral line forms lower in the 



atmosphere (closer to the photosphere), and finally, the op- 
tical radius increases less. As a consequence, the velocity of 
the optical layer decreases for a long-period Cepheid, while 
the gas velocity is supposed to be the same. This picture 
can be generalized at contraction and for the photospheric 

layer. Thus, fo-g = -^^ decreases with the period of the 

Cepheid. Even if the Pfo-g relation seems secured, we keep 
in mind that it must be studied in detail going further in 
our understanding of the dynamical structure of Cepheid's 
atmosphere. Important links between fo-g, the 7- velocity, 
line asymmetry, and velocity gradients should be found. 



5. Direct measurement of velocity gradients from 
HARPS observations 

This section only deals with observations. For the eight 
stars observed, we derive RVc and the line depth as a func- 
tion of the pulsation phase for all spectral lines of Table 
[T] Corresponding uncertainties are estimated based on the 
signal-to-noise ratio. 
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Fig. 4. /S.RVc as a function of the depth D of the spectral line considered. Uncertainties are indicated. Stars are presented 
with increasing period. The 5 Cep model is also indicated for comparison. Linear correlations are derived for all stars 
(see Tab. [3]). The /grad quantity is derived from these relations. 



From interpolated curves (through periodic cubic spline 
functions), we derive Ai?V; = MAX(i?V;)- MIN(i?T4) for 
all stars and lines. We also determine the pulsation phase 
corresponding to the minimum extension of the stars in 
order to derive D. Ai?!^ is represented as a function of 
line depth D for all stars in Fig. [4l Significant linear rela- 
tions (see Table [3]) are found between these two quantities. 
The amplitude of the velocity curves increases with the line 
depth (or with the position of the line- forming region in the 
atmosphere). This validates a posteriori the use of the line 
depth as an estimator of the line-forming region. From the 
coefficients provided in Table[3l we are now able to derive 
the /grad quantity for each star using Eq. [3l 

However, if one wants to compare the velocity gradient 
within the atmosphere of different Cepheids, the question 
of methodology arises. Two strategies are possible. One can 



consider (1) the same line depth for all stars or (2) the 
same spectral line. Because a given spectral line does not 
have the same depth for all stars, the choice is important. 
We first use strategy (1) to determine which spectral lines 
are most fitting for the IBW method, and then consider a 
specific spectral line to compare the velocity gradient in the 
Cepheids' atmosphere. 

First, we test the existence of a linear relation between 
/grad — ^\fj and the logarithm of the period (log P) for 
eleven different line depths from to 0.5 in increments of 
0.05. By a minimization process, we obtain the relations: 
/g,3^d = S{D) \ogP + ZP{D), where the slope {S) and the 
zero-point [ZP] are related to the line depth considered. 
The variables S, ZP, and the reduced x^ are represented 
as a function of the line depth in Fig. [S] The physical in- 
terpretation of these curves is as follows. 
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Fig. 3. fo-g as a function of the period 

Table 3. Linear relations between the amplitude of the 
velocity curves and the line depth (at minimum extension of 
the star), ARVc = a^D + 6o are given for all stars, together 
with the la uncertainty. The reduced x^, defined as Xrod = 

j^_^ with N the number of spectral lines and v the number 
of degrees of freedom is also indicated. 



Star 



ao 



6o 



2 
Arcd 



RTrA 



2.50 ± 4.55 28.48 ± 0.90 



S Cru 


2.13 ± 3.61 


33.33 ± 0.90 


2 


YSgr 


5.76 ± 3.53 


34.01 ± 0.12 


20 


P Dor 


0.86 ± 1.31 


31.59 ± 0.40 


9 


C Gem 


1.09 ± 0.89 


25.35 ± 0.31 


7 


RZ Vel 


8.32 ± 5.95 


47.31 ± 1.02 


6 


£ Car 


2.89 ± 2.26 


32.48 ± 0.67 


15 


RS Pup 


5.89 ± 5.58 


42.19 ± 0.88 


4 



First, for a line depth of zero, the velocity gradient is 
also zero for all stars (or /grad = !)■ Thus, 5 = 0, ZP = 1, 
and there is no uncertainty on /grad- The reduced x^ is 
since the linear relation is perfect. 

Then, for low depths, the /grad estimator is quite dif- 
ferent for each star, but the effect is low and the /grad = 
S{D)\ogP -I- ZP{D) relation is still not very sensitive to 
the velocity gradient. Moreover, the uncertainty on /grad is 
low, and the same is true for S and ZP. The reduced x^ is 
good, as the linearity is well-conserved. 

However, for large depths, the /grad estimator becomes 
more and more sensitive (but also more uncertain) to the 
velocity gradient. The dispersion between the different stars 
is amplified. As a consequence, the linearity is verified at a 
lower level: the reduced x^ gets larger. 

From this picture, our objective is to find the best spec- 
tral lines to use in determining /grad (and then the projec- 
tion factor) in the context of the IBW method. Spectral 
lines with depths lower than 0.1 seem to be the best choice. 
Such spectral lines are indeed less sensitive to the velocity 
gradient. It is obvious that a spectral line that forms close 
to the photosphere implies small differences (in velocity) 
between the line-forming region and the photosphere. In 
this case, /grad is close to 1 and the corresponding uncer- 
tainty is low. Indeed, from Fig. [5l we clearly see that the 
uncertainties on S and ZP are decreasing for low depths. 
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Fig. 5. Slope and zero-point of the Period- /grad relation 
(/grad = S{D)\ogP + ZP{D)) as a function of the line 
depth. The reduced x^ is also indicated. 

From TablcIH we find that the Fel 4896.439 A spectral 
line is the best. It has the lowest depth, averaged over all 
pulsation phases for all stars {D = 0.08). 

In Fig. [6] we present /grad as a function of the logarithm 
of the period for the Fe I 4896.439 A spectral line. We ob- 
tain a linear relation between the velocity gradient and the 
logarithm of the period. Results are given in Table S] for 
observations only (O), and for observations -I- the 5 Cep 
and t Car hydrodynamical models (0-l-C). Consequently, 
we can conclude that the velocity gradient is larger in long- 
period Cepheids than in short-period Cepheids. Moreover, 
if we compare this observational result to the hydrodynami- 
cal models of 5 Cep and (. Car, we find very good agreement 
(see Fig. ^ . The two models presented are thus extremely 
good. 



6. A period-projection factor relation 

We now determine the projection factors of Cepheids con- 
sidering the Fel 4896.439 A spectral line. The Po-factors 
are not determined directly using the continuum intensity 
distribution of the hydrodynamical models. As explained 
in Sect. 14. 4[ these models (except for 6 Cep and £ Car) 
are only used to derive fo-g] no radiative transfer in the 
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Table 5. Derived projections factors for all stars computed from the decomposition presented in Eq. [51 



Name 


HD 


P(b) 
[days] 


Po (c) 


/grad (d) 


fo-3 (C) 


P(f) 


RTrA 


135592 


3.38925 


1.396±o.oio 


0.995±0.009 


0.967±o.oo5 


1.34±o.o3 


S Cru 


112044 


4.68976 


1.392±o.oio 


0.996±o.oo7 


0.966±o.oo5 


1.34±o.o3 


YSgr 


168608 


5.77338 


1.387±o.oio 


0.991±o.oo5 


0.962±o.oo5 


l-32±o.o2 


13 Dor 


37350 


9.84262 


1.380±o.oio 


0.997±o.oo4 


0.955±o.oo5 


1.31±0.02 


C Gem 


52973 


10.14960 


1.380±o.oio 


0.994±o.oo5 


0.953±o.oo5 


1.31±0.02 


RZ Vel 


73502 


20.40020 


1.375±o.oio 


0.994±o.oo4 


0.951±o.oo5 


1.30±o.o2 


£ Car 


84810 


35.55134 


1.366±o.oio 


0.989±o.oo5 


0.944±o.oo5 


1.27±o.o2 


RS Pup 


68860 


41.51500 


1.360±o.oio 


0.995±o.oo5 


0.943±o.oo5 


1.28±o.o2 


S Cep (a) 


213306 


5.419 


1.390±o.oio 


0.990±o.oo5 


0.963±o.oo5 


1.33±o.o2 


£ Car (a) 


84810 


35.60 


1.366±o.oio 


0.988±o.oo5 


0.944±o.oo5 


1.27±o.o2 



5 Cep and £ Car are hydrodynamical models 

The corresponding Julian dates (To) can be found in Paper II. 

Po is derived from the linear limb-darkening laws of Claret et al. (2000) based on the static models of Kurucz (1992). We then 

apply a slight correction based on the S Cep and £ Car hydrodynamical models: po [hydro] = po[geo] — (0.0174 log P — 0.0022) to 

take the dynamical structure of the Cepheid's atmosphere into account. 

/grad is derived directly from observations using Eq|3] It is important to notice that the results indicated here correspond to 

the Fel 4896.439 A line. In the case of a modeled star, it is derived directly from the hydrodynamical model (see Sect. l4.3"1l . 

fo-g is derived directly from the hydrodynamical models (see Sect. l4.4"|l . 

p-factors defined by p = Po/grad/o-g. Po and fo-g are derived from geometrical and hydrodynamical models respectively, /grad 

is derived from observations. 
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Fig. 6. /grad in the case of the Fe I 4896.439 spectral line 
is represented as a fmiction of the period. Open circles cor- 
responds to the hydrodynamic models. The slope and the 
zero-point of the relation depend on the line depth (see 
Fig.©. 

Table 4. Linear relation between /grad and the logarithm 
of the period : /grad = SlogP + ZP for the Fel 4896.439 
A spectral line. 



ZP 



Xi^ed 



O 

o + c 



-0.003 
-0.003 



0.005 
0.004 



0.997 
0.996 



± 0.006 
± 0.005 



atmosphere is computed, and the observed spectral lines 
profiles (and radial velocity curves) are not used to con- 
strain the models. Consequently, we are not certain about 
the dynamical structure of the modeled Cepheids' atmo- 



sphere and, in particular, about /grad- And if /grad is wrong, 
the limb darkening within the line (and thus the po-factor) 
could be affected. Consequently, we consider a conserva- 
tive approach when deriving po [geo] from the hydrodynamic 
parameters of Table [H We then apply a slight correction: 
Po [hydro]=po[geo]-(0. 0174 logP-0. 0022) using the results 
based on d Cep and i Car hydrodynamic models, which 
were studied in detail (Sect. [3]). The resulting linear rela- 
tion is 



Po = [-0.031 ± 0.008] logP + [1.413 ± 0.009]. (5) 

The /o-g quantities have been deduced directly from 
the hydrodynamical models (see Sect. 14. 4|) . and /grad was 
determined from observations in the previous section. From 
Po, /o-g, and /grad, we can now determine consistent pro- 
jection factors for all stars. A summary of the results is 
given in Tableland illustrated in Fig. [71 

By combining all quantities (po, /grad and /o-g), we are 
able to derive a Pp relation for the first time (see Fig. [71- 
b). This result is a combination of observed and theoretical 
considerations. The resulting linear law (including observed 
and modeled stars) is: 



p == [-0.064 ± 0.020] logP -I- [1.376 ± 0.023]. (6) 

This relation holds for the Fel 4896.439 A spectral fine, 
which presents the lowest line depth in our sample. 

From this relation, two facts must be pointed out. First, 
the Ppo relation is based on the general physical properties 
of Cepheids (effective temperature, surface gravity), while 
-P/grad is derived directly from observations and states that 
the velocity gradient within the atmosphere is larger in 
long-period Cepheids. The linearity of the P/o-g relation 
is qualitatively understood (Sect. 14. 3p and has been quan- 
titatively verified based on eight models. Second, results 
concerning the S Cep and i Car hydrodynamic models are 
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Fig. 7. The projection factor {p ~ Po/grad/o-g) as a func- 
tion of the logarithm of the period (diag. b), together with 
Po (diag. a). Black points and open circles correspond to 
observations and models. 



highly secured: (1) the projection factor of S Cep was con- 
firmed observationally by Merand et al. (2005); (2) the ve- 
locity gradient in the atmosphere of these stars is confirmed 
by HARPS observations (see Sect. [5]); (3) the Po estimations 
are coherent at a 0.025 level with the geometrical models; 
and (4) from points 1, 2, and 3, we can reasonably feel se- 
cure in the estimations of /o-g for S Cep, £ Car, and for all 
stars. 



7. Discussion 

The derived Pp relation will be useful in the context of the 
IBW and SB methods. For example, if we compare Eq. [5] 
with the usual value widely used in the community p =1.36 
(Burki et al. 1982), we obtain a correction for the projection 
factor depending on the period. It is then possible to trans- 
late it into a bias on distances and absolute magnitudes. 
By this process, we obtain the relation: 



AMv =G.101ogP-0.03 



(7) 



where AMy is the correction to consider on the PL re- 
lation. We thus conclude that one can make an errors of 
0.10 and 0.03 on the slope and zero-point of the PL rela- 
tion, respectively, if p = 1.36 is used for all stars instead of 
the Pp relation. This correction is, however, only indicative 



because it is indeed restricted to our definition of the pro- 
jection factor (Eq. [T]) and to the Fel 4896.439 A spectral 
line. 

It is now possible to refine the IBW and SB methods. 
First, we suggest using the RVc radial velocity to avoid 
bias related to the rotation velocity of the star (even if 
Cepheids are supposed to be slow rotators) and the width 
of the spectral line. One then has to determine the RVc 
curve and force the average to be zero in order to avoid 
7- velocity effects. Due to our careful definition of p (Eq.[T]), 
the projection factors proposed in this paper are indeed in- 
dependent of the 7-velocity. The spectral line considered 
must have a depth lower than 0.1 and should be the same 
for all considered Cepheids. The low depth of the spectral 
line is required to diminish the impact of the velocity gradi- 
ent. If the Fel 4896.439 A is used, one can use Eq. [6] directly 
to determine the dynamic projection factors of Cepheids. 
If not, we propose the following method. Given the line 
depth of the spectral line considered for each Cepheid, it 
is possible to determine the /grad from Table[3] and Eq. [31 
If the Cepheid being studied is not in our sample. Fig. [S] 
can be used. Then Po can be determined using a geometri- 
cal model. However, the consistency (at a level lower than 
0.025 onp) between interferometric observations, geometri- 
cal, and hydrodynamical models should be studied in detail 
in the future. For /o-g, one can use Eq. 31 even if this re- 
lation has to be confirmed observationally in the future. 
For this purpose, the development of theory and hydrody- 
namical models is required. Finally, the projection factor of 
Cepheids, following our decomposition, is p = Po /grad /o-g- 
This procedure should be applied to avoid bias in the cali- 
bration of the PL relation. 

However, we know that the masking cross-correlation 
method is widely used to increase the signal-to-noise ra- 
tio on radial velocity measurements. In that case however, 
one cannot exclude the impact of the rotation, the spec- 
tral lines' width, and 7-velocities effects. Nevertheless, we 
can still provide a Pp relation that is more appropriate 
considering an average line depth of D = 0.25. We find 
p = [-0.075 ± 0.031] logP -t- [1.366 ± 0.036]. 

Another important point is that we provide visible pro- 
jection factors that should be used with visible spectro- 
scopic observations. If one used infrared spectroscopic ob- 
servations to derive the pulsation velocity, one should use 
specific infrared projection factors. Indeed, in the infrared, 
the limb darkening is supposed to be lower and the corre- 
sponding po-factors higher (certainly about 4%). But, spec- 
tral lines also form higher in the atmosphere (i.e. in the up- 
per part of the atmosphere), which supposes a lower /grad- 
More studies have to be carried out to derive an infrared 
Pp relation. 



8. Conclusion 

In the application of the IBW method, the projection fac- 
tor is a key quantity. Up to now, the period-dependency 
of the projection factor has never been studied in detail. 
Here, we have presented a new spectroscopic method for 
directly measuring the velocity gradient in the Cepheids' 
atmosphere. This method has been successfully validated 
by the hydrodynamical models of S Cep and i Car. We find 
a physical relation between the period of the star and its 
dynamical atmospheric structure. 
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The models also show that the optical layers (observed 
by continuum interferometer) and the gas layers have to be 
distinguished in the interferometric definition of the pro- 
jection factor. However, this quantity is still very difficult 
to determine directly from observations. 

Combining the results obtained directly from observa- 
tions and our knowledge of the dynamical structure of the 
S Cep and £ Car atmosphere, we have been able to derive a 
very consistent Pp relation for the Fel 4896.439 A spectral 
line: 

p = [-0.064 ± 0.020] logP + [1.376 ± 0.023]. (8) 

We emphasize that, if a constant projection factor is 
used to constrain the PL relation, an error of 0.10 and 0.03 
magnitudes can be done, respectively, on the slope and zero- 
point of the PL relation. This can even be much more if 
the wrong definition of the radial velocity is used or if one 
does not consider 7-velocity effects. We have thus presented 
(see discussion) a careful methodology to be applied in the 
context of the IBW and SB methods. 
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